Specific heat anomalies of open quantum systems 



Gert-Ludwig Ingold, 1,2 ^ Peter Hanggi, 1 and Peter Talkner 1 

1 Institut fur Physik, Universitat Augsburg, D-86135 Augsburg 
2 Laboratoire Kastler Brossel, CNRS, ENS, UPMC, 
Campus Jussieu Case 74, F-75252 Paris Cedex 05, France 

The evaluation of the specific heat of an open, damped quantum system is a subtle issue. One 
possible route is based on the thermodynamic partition function which is the ratio of the partition 
functions of system plus bath and of the bath alone. For the free damped particle it has been shown, 
however, that the ensuing specific heat may become negative for appropriately chosen environments. 
Being an open system this quantity then naturally must be interpreted as the change of the specific 
heat obtained as the difference between the specific heat of the heat bath coupled to the system 
degrees of freedom and the specific heat of the bath alone. While this difference may become 
negative, the involved specific heats themselves are always positive; thus, the known thermodynamic 
stability criteria are perfectly guaranteed. For a damped quantum harmonic oscillator, instead of 
negative values, under appropriate conditions one can observe a dip in the difference of specific heats 
as a function of temperature. Stylized minimal models containing a single oscillator heat bath are 
employed to elucidate the occurrence of the anomalous temperature dependence of the corresponding 
specific heat values. Moreover, we comment on the consequences for the interpretation of the density 
of states based on the thermal partition function. 



I. INTRODUCTION 

An open classical system in contact with a heat bath 
can often be modeled in terms of a Langevin dynamics 
with constant friction and white Gaussian noise sources 
obeying a fluctuation-dissipation theorem [TJ [3J |3J . A re- 
markable feature then is the circumstance that the equi- 
librium statistics of the open classical system turns out 
to be independent of the coupling strength between the 
system and the heat bath. In other words, the canonical 
equilibrium for a classical damped Langevin dynamics 
agrees with the canonical equilibrium of the isolated sys- 
tem. This feature is rooted in the fact that in this case the 
so termed "Hamiltonian of mean force" is still given by 
the bare system Hamiltonian [4]. In clear contrast, this 
property in general no longer holds true for open systems 
in the quantum regime beyond the weak-coupling limit 
[5]. In particular, the canonical equilibrium state of an 
open quantum system then typically involves an explicit 
dependence on the system-bath coupling strength. 

Motivated by this fact, the study of the specific heat 
beyond the weak-coupling limit has recently received con- 
siderable attention, in particular in view of the validity of 
the Third Law of thermodynamics [BJ HI EO EU EH EH • 
Apart from fundamental thermodynamical questions the 
study of the specific heat in the quantum regime is also of 
interest because it can be related to entanglement prop- 
erties [IB] , 

Recently, two different routes towards the evaluation 
of a specific heat were proposed and discussed [H H] . One 
possibility is based on the thermal expectation value of 
the Hamiltonian describing the isolated system. Another 
approach, on which we will focus in this paper, starts 



'Electronic address: gert.ingold@physik.uni-augsburg.de 



out from the thermodynamic partition function of the 
dissipative system [6l[g[l4l[I32[3Il[I3[ig 



Z = 



Tr s+B [exp(-/3g)] 
Tr B [exp(-/3ff B )] 
where the total Hamiltonian 



H — Hq + Hvt + H\ 



SB 



(1) 



(2) 



consists of terms describing the system, the bath, and 
the system-bath coupling, respectively. In the absence of 
a coupling between system and bath, Z reduces to the 
partition function of the system. The partition function 
([lj appears naturally in the Feynman- Vernon approach 
to dissipative systems [15j EBJ Ell EHj and can be related 
to equilibrium properties of the system [XT] . 

From 0, one obtains by means of standard thermo- 
dynamic relations a specific heat [25], reading 



C=k B P 2 ^HZ). 



(3) 



Here, fc B is the Boltzmann constant and the temperature 
T appears through /? = l/fc B T. 

In the following, we will assume the bath to consist of 
harmonic oscillators and the coupling to be bilinear in 
system and bath coordinates 5 . In such a framework 
Ref. [5] found for the damped free particle, that under 
certain circumstances the specific heat ([3| can become 
negative. In the case of the Drude model, for example, 
where the Laplace transform of the damping kernel is 
given by 



7(2) = 



7^d 
z + ud 



(4) 



the specific heat exhibits negative values at low temper- 
atures if the damping constant 7 exceeds the cut-off fre- 
quency wd, i.e. 7 > wd- This behavior is depicted in 
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FIG. 1: The specific heat as defined inj3| for a free damped 
particle with a Drude damping kernel (Ml is shown as a func- 
tion of temperature for cuo/'y = 0.2, 1, and 5. Note that in 
contrast to Ref. J5J , temperature is given in units of the cutoff 
frequency uu instead of the damping strength 7. 



Fig. [T] For a general damping kernel, a negative specific 
heat will appear at low temperatures if 7'(0) < —1, see 
Ref. [9]. No negative specific heat was found if the free 
particle is replaced by a harmonic oscillator. 

In this paper, we will elucidate the origin of the nega- 
tive specific heat for the case of the free damped particle 
and its absence for the damped harmonic oscillator by 
considering two minimal models. Before doing so, a gen- 
eral comment on the appearance of a negative specific 
heat is appropriate. 

The rationale behind the definition of a specific heat 
based on the thermodynamic partition function ([TJ is 
that this partition function should be associated with an 
effective description of an open system. From this point 
of view, a negative specific heat may appear disturbing 
because it raises doubts as to the thermodynamic stabil- 
ity of the system. 

However, the meaning of the thermodynamic partition 
function ([lj of the open system can be better understood 
from the point of view of its finite coupling to the heat 
bath. The specific heat ^ can in fact be expressed as 



C — Cs+b — C) 



(5) 



where Cs+b is the specific heat of system and heat bath 
while Cb is the specific heat of the heat bath alone. 

Therefore, C describes the change of the specific heat 
when the heat bath is enlarged by coupling it to sys- 
tem degrees of freedom. This difference expression in Q 
must be expected on physical grounds when dealing with 
an open system that is not heat-isolated from its environ- 
ment. For example, take a certain amount of an agent 
within a container: The established experimental proce- 
dure to determine the specific heat of this agent is first 
to measure the specific heat of the empty container and 
to subtract this value from the measured specific heat 
of the combined system to finally arrive at the specific 



heat of the agent alone. For a macroscopic amount of the 
agent this value is truly agent specific, i.e. independent of 
the particular nature of the container and its interaction 
with the agent (provided the interaction of the agent with 
the container is short ranged). For a nanoscopic system, 
however, the energy that is contained in the system-bath 
coupling typically cannot be neglected and consequently 
will influence the specific heat of such a system. Notably 
the thermodynamic internal energy of the open system 
also typically differs from the thermal expectation of the 
bare system Hamiltonian [HE]. As a consequence, the 
values of the specific heat evaluated along such different 
routes then differ as well [HE]- 

The coupling can thus result in a negative specific heat 
of the open quantum system. However, the involved spe- 
cific heats, i.e. those of system plus bath on the one hand 
and of the bath alone on the other hand, are each positive 
so that no issues concerning the thermodynamic stabil- 
ity arise |26] . Such physical situations are not uncommon 
as evidenced by recent discussions in the context of the 
Casimir effect [27], the multichannel Kondo effect (28], or 
the physics of mesoscopic superconductors that contain 
magnetic impurities |29j . 



II. FREE PARTICLE 

In order to elucidate the appearance of a negative spe- 
cific heat ([3]) it is sufficient to consider a stylized, minimal 
model where the "bath" consists of only a single degree 
of freedom described by the Hamiltonian 



/b 2 



P 



(6) 



where /b denotes the spring constant. In the following, 
we will study a system governed by the Hamiltonian 



2M 



2 



(7) 



both in the cases of a free particle (spring constant /g = 
0) and of a harmonic oscillator (/g > 0). The coupling 
Hamiltonian is given by 



Hi 



SB 



-fmQ 



h 



Q 



(8) 



where the last term renormalizes the potential in order 
to ensure translational invariance in the case of the free 
particle. Figure [2] illustrates our two minimal models. 
The starting point is a single bath oscillator with mass 
m as depicted in Fig. [2^. Coupling the system mass M 
to the bath oscillator leads to the harmonically coupled 
system of two masses shown in Fig. [2]d. If the system 
degree of freedom corresponds to a harmonic oscillator, 
we obtain the mechanical system of Fig. [2]:. 

It is quite obvious that an environment consisting of 
one degree of freedom does not suffice to replace any re- 
alistic heat bath. In particular, it does not lead to a 
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FIG. 2: (a) One single bath oscillator represented by a mass 
m harmonically coupled to a wall of infinite mass, (b) Bath 
oscillator coupled to a free system degree of freedom, (c) 
Bath oscillator coupled to a harmonically bound system de- 
gree of freedom. The spring constants of the bath and system 
oscillator are denoted by /b and fs, respectively. 



truly dissipative behavior of the open system to which it 
couples. Nevertheless, it turns out that even such min- 
imal bath models give rise to the same thermodynamic 
anomalies that are also encountered with more realistic, 
large environments. 

We begin our analysis with a free particle in contact 
with the singlc-dcgrcc-of-frecdom environment described 
by eqs. Q and System and bath are assumed to 
stay in thermal equilibrium with each other at the inverse 
temperature f3. Hence, the density matrix of the total 
system is given by a Gibbs state reading 



PSB = Z S b exp [-p(H a + H B + H SB ) 



(9) 



where Z$b = Trexp [— (3(Hs + Hb + -Hsb)] denotes the 
partition function of the total system. 

The partition function Zb = Tr exp [— (3Hb] of the iso- 
lated bath degree of freedom is given by 



Zb = 



1 



2 sinh 



where 



1/2 



(10) 



(11) 



is the frequency of the bath oscillator. From (D the spe- 
cific heat of this bath follows as 



C B = k B g\-^- 



(12) 



with the abbreviation 



sinh (a;) 



(13) 



If we add the system degree of freedom in order to ob- 
tain the mechanical system shown in Fig. [2j>, the par- 
tition function contains contributions of two degrees of 
freedom related to the center-of-mass and the relative 
motion. The first one is described by a free particle with 
an effective mass m + M while the second degree of free- 
dom corresponds to a harmonic oscillator with effective 
mass mM/(m + M) and the frequency 



1 



AI 



1/2 



(14) 



As discussed in the introduction, negative values of the 
specific heat occur for the Drude model if the damping 
strength exceeds the cut-off frequency of the heat bath. 
Within our minimal model the specific heat may become 
negative if the mass ratio m/M exceeds a value slightly 
above 4. Then, Q is significantly larger than to, a fact 
which will be relevant for the discussion of the specific 
heat ( 16 1 below. 



In order to obtain a well-defined partition function for 
the free particle, we restrict its motion to a region of 
length L. This length is supposed to be sufficiently large 
such that the energy level spacing can be neglected if 
compared with the thermal energy fc^T [9]. Under this 
condition L will turn out to be irrelevant in the sequel. 

The partition function of system plus bath consists of 
a product of contributions arising from the two normal 
modes, i.e. the center-of-mass and relative motion, and 
thus reads 



^SB 



L 

h 



2n{m + M) 



1/2 



2 sinh 



(hf3u) 

V~2~ 



(15) 



From (10 1, (12 1 ( 13 1 and (|15[ ) it is straightforward to 
evaluate the specific heat (13} which becomes 



C 
k B 



H(3lu 



(16) 



The first term arises from a free particle while the second 
and third term describe the change in specific heat due 
to the increase in the oscillator frequency from uj to uj as 
given by ( 14 1. 

In Fig. 3pthe contributions to (16) are sketched. The 
upper dashed curve corresponds to the first two contribu- 
tions arising from system and bath. It contains the fce/2 
from the isolated free particle and a contribution from 
the harmonic oscillator which is strongly suppressed at 
low temperatures and reaches k B for high temperatures. 
The lower dashed curve corresponds to the third term 
in (16 1. The main point to note is the relative shift in 



temperature of the two contributions due to the change 
of the oscillator frequency. In the presence of the system, 



4 
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of system plus bath are readily obtained as 



FIG. 3: The difference of specific heats as function of the 
temperature for the case of a free particle coupled to a single 
oscillator bath for the mass ratio m/M = fO. The upper 
dashed curve corresponds to the first two terms in \1Q\ , i.e. 
to the specific heat of the system shown in Fig. [2V> while the 



lower dashed curve corresponds to the third term in (16 1, i.e 



to the negative of the specific heat of the bath oscillator shown 
in Fig. [2^. The resulting difference is depicted as solid line 
and displays a temperature region where it takes on negative 
values. 



the oscillator frequency is increased according to (14 1. As 



a consequence, there is a temperature window, where the 
specific heat of system and bath is already significantly 
suppressed while this is not yet the case for the bath os- 
cillator alone. In this regime, the difference (16 1 of the 



specific heats can become negative. Note that in contrast 
to Fig.[T]this temperature window does not extend all the 
way down to zero temperature. This is explained by the 
fact that the bath consists of only one oscillator so that 
the low-frequency oscillators present in the Drude model 
Q on which Fig. [I] is based are missing. 



III. HARMONIC OSCILLATOR 

After having explained the origin of a negative specific 
heat for a free particle in contact with an environment we 
will now address the question why such a scenario cannot 
be realized for a harmonic oscillator. 



A. Minimal environment 

To this end, we determine the difference of specific 
heats for the mechanical systems shown in Figs. [2^i and 
Ec. The system frequency associated with the spring con- 
stant /s is given by Q = (/s/M) 1 / 2 . The eigenfrequencies 



1 /m + M 



-a; 2 + n< 



± 



71/ 

1 / m + M 



M 



uj 2 + n 2 



.1/2 (17) 



> 2 n 2 



One can show that w_ < lu < lj + where w_ = u> or 
lu + = uj for m — depending on whether cu/fl > 1 
or uj/CI < 1. For increasing mass ratio m/M and fixed 
frequency ratio u/fl both gaps from u> to w_ and w + 
widen. In the limit m 3> M, at any fixed frequency ratio 
oj/n one finds w_ = (M/m) 1 / 2 ^ and lu + = (to/M) 1 / 2 ^. 
From ([TJ and ^ one then obtains 



C fh(3uj+\ fh(3oj- 
^= 9 {—) +g { — 



fJ 



h(3u 



(18) 



where k-ag is the specific heat of a harmonic oscillator 
defined in Eqs. (12 1 and (13 1. A typical scenario for the 
case of sufficiently well separated frequencies is sketched 
in Fig. [4] Although for the harmonic oscillator a dip in 
the specific heat may appear, no negative values can be 
obtained. The main difference to the case of a free par- 
ticle is the specific heat of the isolated system degree of 
freedom: While for not too low temperatures the specific 
heat for the harmonic oscillator equals /cb it is only half 
as large for the free particle. The difference of the specific 
heats of the bath oscillator in the presence and absence 
of the system degree of freedom may reach values up to 
&b, thereby opening up the possibility of negative values 
of the specific heat for the free particle but not for the 
harmonic oscillator. 



B. Drude bath 

The dip in the specific heat of the harmonic oscilla- 
tor in contact with the minimal environment can also be 
observed for a bath giving rise to ohmic damping, i.e. 
where 7(0) > 0. For the ratio of partition functions 
one finds 



1 00 2 



(19) 



where v n — 2iTn/hp are the Matsubara frequencies. 

We specifically consider a Drude model characterized 
by the damping kernel Q with a cutoff frequency wd- 
Introducing the quantities 



and 



A; 



At 



h(3\i 



2?r 



1,2,3 



(20) 



(21) 
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FIG. 4: The difference of specific heats as function of the 
temperature for the case of a harmonic oscillator coupled to 
a single oscillator bath. The upper dashed curve corresponds 
to the first two terms in flTSp, i.e. to the specific heat of 
the system shown in Fig. I2h while the lower dashed curve 
corresponds to the third term in (TlSp , i.e. to the negative of 
the specific heat of the bath oscillator shown in Fig. [2^,. The 
resulting difference is depicted as solid line and displays a dip 
for the chosen parameter values £1 = u and m/A'I = 10. 



with 

x 3 + ld-qx 2 + (■yuj-D + Q 2 )x + udO 2 
= (x - Ai)(.t - X 2 )(x - A 3 ) , 

we obtain from ([3| with (19 1 

£ = 1 + A?V'(1 - A x ) + A^'(l - A 2 ) 
+ A^'(1-A 3 )-A^'(1-A D ) 



(22) 



(23) 



Here, ip'(x) denotes the trigamma function. The specific 
heat (23 1 is shown in Fig. [5] for 7 = 50 and wd = O.lfi 
where is the frequency of the system oscillator. In con- 
trast to the undamped case, the specific heat increases 
linearly with temperature with a slope proportional to 
the damping constant. Then, for sufficiently small cutoff 
frequency the specific heat goes through a dip at low tem- 
peratures before it asymptotically approaches its high- 
temperature value Kb- 



IV. DENSITY OF STATES 

One can use the partition function to formally de- 
fine an effective density of states p(E) of the system by 
means of the relation I2T1 



Z = 



dEp(E) exp(-f3E) . 



(24) 



For the free particle, it was found in Ref. [9] that the 
density of states obtained from (24 1 can become nega- 



FIG. 5: Specific heat |23f of a damped harmonic oscillator 
for Drude damping with 7 = 5Q and uio ~ 0.1Q. The inset 
shows a close-up of the low-temperature regime. 



model with ujd < 7- For an environment consisting of 
a single bath mode like in the minimal models discussed 
in Sect. [IT]and |III| it was demonstrated, that the appear- 
ance of negative contributions to the effective density of 
states is generic. In particular, the minimal model for 
a harmonic oscillator coupled to an environmental mode 
displayed in Fig. [2J; leads to negative delta functions in 
the density of states, see Ref. |9]- This is in contrast 
to our finding for the specific heat which remains always 
positive for the harmonic oscillator. Interestingly, for a 
heat bath exhibiting a continuous distribution of bath 
oscillators, one obtains again a positive density of states 
despite of the fact that dips in the specific heat like the 
one shown in Fig. [5] can be observed. Fig. [6] displays the 
density of states for the same parameters as employed in 
Fig. [5] 

An important difference between the specific heat and 
the density of states as defined by (24 1 lies in the fact that 



the latter cannot be interpreted in terms of a difference 
of two densities of states. We recall that such an inter- 
pretation was possible for the specific heat only because 
according to ^ it depends linearly on the logarithm of 
the partition function ([lj. The specific heat shares this 
property with other thermodynamic quantities like the 
internal energy and the free energy. The absence of a 



logarithm of the partition function in (24 1 indicates that 



tive for appropriately chosen environments, e.g. a Drude 



the effective density of states does not lend itself to an in- 
terpretation in terms of the difference of two densities of 
states. Despite the fact that the effective density of states 
of a damped harmonic oscillator in the weak-coupling 
limit, i.e. when the damping strength represents the 
smallest frequency scale, displays resonances at the ex- 
pected energies and even yields the correct level widths 
[23], it therefore remains unclear whether the meaning 
of the effective density of states goes beyond that of a 
merely formal notion. 
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FIG. 6: Density of states (24 1 of a damped harmonic oscil- 
lator for Drude damping with 7 = 5Q and ujo = 0.1Q. For 
these parameters, (221 yields a resonance frequency of about 
1.2240. A delta function contribution of weight one present 
at the ground state energy Eo of the damped harmonic oscil- 
lator is not shown [30] . 



a damped harmonic quantum oscillator, the difference 
of specific heats cannot become negative but may dis- 
play a dip instead. The difference in the behavior of the 
free particle and the harmonic oscillator can be traced 
back to the specific heat of the uncoupled system which 
for the free particle is smaller by a factor of two. All 
those quantities that are obtained from the logarithm of 
the partition function by means of a linear operation can 
be interpreted as differences between the corresponding 
quantities of the total system plus bath and of the bath 
alone. This reasoning though does not apply to the den- 
sity of states which is the inverse Laplace transform of 
the partition function itself, see ( 24 ) . By this transforma- 



tion the ratio of two partition functions yields a complex 
quantity that obviously cannot be interpreted as a differ- 
ence of the densities of states of the total system and of 
the environment. 
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